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Abstract 

We exhibit infinite families of planar graphs with real chromatic 
roots arbitrarily close to 4, thus resolving a long-standing conjecture 
in the affirmative. 

1 Introduction 

The chromatic polynomial P(G, x) of a graph G is the function that counts 
the number of proper x-colourings of a graph when x is a positive integer. 
It was introduced by Birkhoff [6] in 1912, in the hope that a quantitative 
study of the numbers of colourings of planar graphs could be used to show 
that P(G, 4) > 0, thus resolving the 4-colour conjecture. Although defined 
on positive integers, it is well-known that P(G,x) is a polynomial in x, 
and therefore we can view it as a function on the real or complex numbers. 
Although they were unable to prove the 4-colour conjecture with chromatic 
polynomials, Birkhoff and Lewis [5] showed that for G planar, P(G,x) > 
for all x G [5, oo), and they conjectured that P(G,x) > for all x G [4, oo), 
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thus providing the first results and problems in the theory of real chromatic 
roots. 

With the 4-colour conjecture still unsolved, the study of the complex chro- 
matic roots of planar graphs was initiated in the 1960s (Hall, Siry & Vander- 
slice [10], Berman & Tutte [3]) when the availability of computers made em- 
pirical study of complex chromatic roots feasible. A decade later, at around 
the same time as the 4-colour conjecture became a theorem with a decidedly 
combinatorial proof, Beraha & Kahane [2] found families of planar graphs 
with complex chromatic roots arbitrarily close to 4, thus providing strong 
evidence that there could never be a complex-analytic proof of the 4-colour 
theorem. 

Despite this, the study of real and complex chromatic roots of not- necessarily- 
planar graphs has developed into a substantial theory, significantly enhanced 
by the relationship between the chromatic polynomial and the partition func- 
tion of the g-state Potts model of statistical mechanics. Jackson [11] provides 
an excellent and comprehensive survey of the main results and questions in 
the theory of chromatic roots, while an accessible introduction to the connec- 
tions between the chromatic polynomial and the Potts model can be found 
in Sokal [17]. 

One of the unresolved problems described by Jackson [11] was a conjecture of 
Beraha which (among other things) asserted the existence of planar graphs 
with real chromatic roots arbitrarily close to 4. At the time of Jackson's 
survey the largest known real chromatic root of a planar graph was 3.8267 . . . 
obtained from a 21-vertex graph found by Woodall. In this paper we resolve 
this conjecture in the affirmative, describing a number of related families of 
planar triangulations with real chromatic roots arbitrarily close to 4. 

2 Double-ended lattice graphs 

Some of the most important families of graphs from a statistical mechanics 
point of view are rectangular subsets of infinite two-dimensional lattices, 
particularly the square and triangular lattices (see Figure 1 for the triangular 
lattice). 

Physicists are usually interested in the limiting behaviour of the complex 
chromatic roots as the strip becomes infinitely wide and tall. In practice, this 
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Figure 1: A 4 x 6 strip of the infinite triangular lattice 



is studied by fixing one of the parameters, say the width, and examining the 
curves where the chromatic roots accumulate as the height tends to infinity. 
Repeating this calculation for a number of different fixed widths then gives 
insight into how the entire curves move as the width increases. In addition, 
they consider the cases of periodic boundary conditions where the left edge of 
the strip is viewed as being adjacent to the right edge as though the strip were 
wrapped around a cylinder, and free boundary conditions otherwise. These 
different boundary conditions are denoted by using P and F respectively 
as subscripts, so for example mp x np denotes an m x n rectangular strip 
(of whichever type of lattice is under consideration) with periodic transverse 
(i.e. left-right) boundary conditions and free longitudinal (i.e. top-bottom) 
boundary conditions. 

For families of graphs of this form with one parameter fixed, the limiting 
behaviour of the chromatic roots can be explained by a theorem of Beraha, 
Kahane & Weiss [1] which shows that the accumulation points of the roots 
form one or more curves and isolated points in the complex plane. The lo- 
cation and nature of these curves has been heavily studied by a number of 
researchers, in particular the prolific group of Shrock, Chang & Tsai who, 
with various collaborators, have studied both the full Tutte polynomial (a.k.a 



"the partition function of the g-state Potts model" ) and the chromatic poly- 
nomial for various different lattices with a number of different boundary 
conditions. Shrock [16] gives a good introduction to their earlier work on 
chromatic polynomials, while the recent work [8] contains many other refer- 
ences. For the triangular lattice, which is the focus of this paper, the full 
Tutte polynomial for strips of this lattice has been computed and studied for 
widths up to 5 by Chang, Jacobsen, Salas & Shrock [7], and the chromatic 
polynomial for widths up to 10 by Jacobsen, Salas & Sokal [12]. 

The graphs found by Beraha & Kahane [2] with complex chromatic roots 
arbitrarily close to 4 are 4p x np strips of the triangular lattice with an ad- 
ditional vertex joined to each end. Rocek, Shrock & Tsai [13, 14] considered 
the effect of adding a number of different end graphs to either one end or 
both ends of a strip, and it is these double-ended lattice graphs that we can 
use to find real chromatic roots arbitrarily close to 4. 

In particular we show that if we attach suitable pairs of planar end-graphs 
to the top and bottom of the cylindrical triangular lattice 4p x uf, then the 
resulting families of graphs have real chromatic roots arbitrarily close to 4. 

3 Chromatic Polynomials 

Let A and B be two graphs each with a distinguished 4-cycle, say 01020304 
and 61 02&3&4 respectively, and suppose that G is the graph obtained by glu- 
ing together A and B, i.e. identifying Oj with fej for 1 < % < 4. Any proper 
colouring of G determines proper colourings of A and B such that the dis- 
tinguished 4-cycle is coloured identically in both graphs. Conversely, if we 
have enough information about the colourings of A and B and the colourings 
of the distinguished 4-cycle in both A and B, then we can determine the 
chromatic polynomial of G. 

Given a graph A with distinguished 4-cycle 01020304, we can assign one of 
four types to each proper colouring if of A according to the colour partition 
it induces on the 4-cycle, as follows: 

Type 1 </?(oi) = (p(a 3 ) and ip(a 2 ) = (p(a 4 ), 
Type 2 ip(a 1 ) = </?(o 3 ) and ip(a 2 ) ^ y?(a 4 ), 
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Type 3 </?(ai) ^ ip(a 3 ) and ip(a 2 ) = <p(a 4 ), 



Type 4 ip^) ^ </?(a 3 ) and </?(a 2 ) 7^ y?(a 4 ). 



The number of colourings of each type is just the chromatic polynomial of 
an auxiliary graph obtained from A by identifying two vertices if they are 
required to have the same colour, and adding an edge between them if they 
are required to have different colours. Let Pi(A, x), P 2 (A,x), P 3 (A,x) and 
P^(A,x) denote the four chromatic polynomials corresponding to the four 
classes above, and note that 

P(A, x) = Pi(A, x) + P 2 (A, x) + P 3 (A, x) + P 4 (A, x). 

We will find it convenient to express this as a vector in Z[x} 4 and so we define 
the partitioned chromatic polynomial Q(A,x) of graph with a distinguished 
4-cycle to be 

/ Pi(Ax) \ 

P2(A,x) 
P 3 (A,x) 
V P4(A,x) J 



Q(A,x) 



3.1 Lemma. Let A and B he two graphs with distinguished 4-cycles and 
corresponding partitioned chromatic polynomials Q(A,x) and Q(B,x). Then 
the chromatic polynomial of the graph G obtained by gluing together A and 
B is the sole entry of the lxl matrix 



Q(A) T DQ(B) 



D 



( l/{xh 



\ 





l/(x) 3 

l/(x) 3 





\ 



where 





1/(Z> 4 / 

and {x)k denotes the k'th falling factorial x(x — 1) • • • (x — k + 1). 



Proof. For a positive integer x, the product Pi(A, x)Pi(B, x) counts the 
number of pairs (p, ip) such that p is a proper re-colouring of type i of A, 
and ip is a proper re-colouring of type i of B. The two colourings can be 
combined to form a colouring of G if and only if they both use the same 
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Figure 2: L is one layer of the cylindrical triangular lattice 



colours on the distinguished 4-cycle, and all proper x-colourings of G arise 
in this fashion. If the colourings of type % use s colours on the distinguished 
4-cycle, then the fraction of these pairs of colourings using the same colours 
on the distinguished 4-cycle is exactly l/(x) s . Summing the four terms of 
this form, and expressing the result in matrix terms gives the stated result, 
n 

Now consider the "gadget" L shown in Figure 2; this is actually a 4px2^ strip 
of the triangular lattice, or in other words a single "layer" of the cylindrical 
triangular lattice of width 4. Let M denote the matrix whose rows and 
columns are indexed by the 4 types of colour partition of a 4-cycle and where 
the entry M^- is the polynomial counting the number of x-colourings of L that 
are of type % on the outer 4-cycle, and type j on the inner 4-cycle. Every 
partition of L into s independent sets contributes a term of (x) s to one of 
the entries of M, and a simple computation yields 

/ {x)a (x)5 (x)5 (x)g \ 

M= {x) 5 (x) 4 + 2(x) 5 + (x) 6 {x) 4 + 2{x) 5 + {x) 6 4{x) 5 + A{x) 6 + {x) 7 
(x) 5 (x> 4 + 2(x> 5 + (x) 6 (x> 4 + 2{x) 5 + (x) 6 4(x) 5 + 4(x) 6 + (x) 7 
V (x) 6 4(x) 5 + 4(x) 6 + (x) 7 4(x) 5 + 4(x) 6 + (x) 7 M 44 / 

where 

M 44 = 2(x) 4 + 16(x) 5 + 20(x) 6 + 8(x) 7 + (x) 8 . 
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Now we consider the effect of adding a layer of the cylindrical triangular 
lattice to an existing graph. In particular, take a graph A with a distinguished 
4-cycle, and identify this 4-cycle with the inner 4-cycle of L, forming a new 
graph A' whose distinguished 4-cycle is the outer 4-cycle of L. In order to 
repeat this operation with A', we need some way to "transfer" the colouring 
information from the inner 4-cycle to the outer 4-cycle. 

3.2 Lemma. Let A be a graph with a distinguished 4-cycle, and suppose 
that A' is obtained by identifying the inner 4-cycle of L with the distinguished 
4-cycle of A, and declaring the outer 4-cycle of L to be the distinguished 4- 
cycle of A'. Then the partitioned chromatic polynomial Q(A',x) is given 
by 

Q(A') = MDQ(A). 

Proof. This is a counting argument analogous to that of Lemma 3.1, where 
the proper colourings of A' are counted according to the type of colour par- 
tition on the inner 4-cycle of L, and then allocated to the result according 
to the type of colour partition induced on the outer 4-cycle of L. □ 

The matrix MD is the transfer matrix for the cylindrical triangular lattice, 
because it transfers the colouring information from one layer of the lattice 
to the next. In one form or another transfer-matrix techniques are the main 
tool for studying chromatic and Tutte polynomials of lattice graphs — more 
details on two approaches to using transfer matrices can be found in Chang 
k Shrock [8] and Salas & Sokal [15]. 

3.3 Theorem. Let A and B be graphs with distinguished 4-cycles, and let 
XA,B(n) be the graph obtained from the triangular lattice strip 4p x uf by 
gluing A to the 4-cycle at the top of the strip and B to the 4-cycle at the 
bottom. Then the chromatic polynomial of XA,B{n) is the sole entry of the 
lxl matrix 

Q(A) T D(MD) n - l Q(B). 
Proof. This follows directly from Lemma 3.1 and Lemma 3.2. □ 

Finally, we observe that Woodall's graph with real chromatic root 3.8267 . . . 
is the graph X HjWi (2) where H and W 4 = C 4 + K x are the graphs shown in 
Figure 3 and Figure 4 respectively with the outer face being the distinguished 
4-cycle in each case. 
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Figure 3: H is a triangulation of a square 




Figure 4: The 4-wheel W 4 = C 4 + fTi 
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1 


538 
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—28 


—5244 


—5244 


—5246 


X 7 


362 


35078 


35078 


35148 


X 6 


-2846 


-169490 


-169490 


-170548 


X 5 


15036 


604806 


604806 


613920 


x A 


-55448 


-1595807 


-1595807 


-1645010 


X 3 


142716 


3051803 


3051803 


3222645 


x 2 


-246724 


-4024676 


-4024676 


-4397342 


x 1 


258889 


3286881 


3286881 


3753346 


1 


-124884 


-1255163 


-1255163 


-1511254 



Table 1: Components of the partitioned chromatic polynomial of H . 



4 The family Xh,wJ<p) 



In this section we give computational results on the real roots of the family 
Xh, W4,( n ), which although not quite the smallest family having real chromatic 
roots arbitrarily close to 4, is the first family for which this was noticed. It 
is immediate that 



Q(w 4 ,x) 



( (x) s \ 
(x) 4 
(x) 4 



and Q{H) is shown in Table 1, where each term is given as the product of a 
simple multiplier and an explicitly listed polynomial of degree at most 11. 

Using the expression given in Theorem 3.3 and a suitable computer algebra 
package, it takes only minutes to calculate the chromatic polynomials for 
-^if,w 4 (l) to X HjWi (100). The roots of these polynomials were then computed 
and the largest real roots listed in Table 2. 

For these polynomials, it is also possible to compute their complex chromatic 
roots; the roots for X H; vk 4 (100) are plotted in Figure 5 and are in close 
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n 


Max. real root tor X h, w 4 \ n ) 


1 


3.7924699360 


2 


3.8267852044 


3 


3.8483432574 


4 


3.8637744449 


5 


3.8756040984 


10 


3.9100811222 


on 

20 


3.9388450668 


30 


3.9524053758 


40 


3.9605533349 


50 


3.9660736920 


60 


3.9700969536 


70 


3.9731775927 


80 


3.9756223799 


90 


3.9776159809 


100 


3.9792767496 



Table 2: Largest real chromatic roots of X H>W4 (n). 

agreement with the theoretical limiting curves as first shown in Beraha & 
Kahane [2]. 

The computation of the complex chromatic roots starts to become time con- 
suming for larger chromatic polynomials, but we can go much further if we 
are only interested a single real root. By repeated squaring, we can easily 
compute (MD) 2 , (MD) 4 , . . ., (MD) 512 and the bisection method (with ra- 
tional arithmetic) can be used to locate a real chromatic root near 4, which 
is presumably the largest. The values found in this manner are given in 
Table 3. 

4.1 Notes on the computations 

As the results of this section are computational in nature, it is important 
to verify them as much as feasible. To this end, the partitioned chromatic 
polynomial of the base graph H was computed by two independent programs, 
being Gary Haggard's fast chromial program [9] and a slower but much sim- 
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Figure 5: Complex chromatic roots of w 4 (100) 
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n 


Value 


2 


3.848343257 


4 


3.875604098 


8 


3.905148525 


16 


3.932717391 


32 


3.955237394 


64 


3.971732100 


128 


3.982848013 


256 


3.989898687 


512 


3.994181944 



Table 3: Real chromatic roots of X HjWi (n + 1) close to 4. 

pier program of the author. The computation of the chromatic polynomials 
was performed symbolically by Maple and then using exact arithmetic we 
confirmed that the polynomial satisfied Tutte's identity for planar triangu- 
lations, i.e that 

P(G, r + 2) = (r + 2)t 3 "- 10 P(G, t + l) 2 

where r is the golden ratio (1 + \/5)/2 and n is the number of vertices of G. 
The computation of the complex chromatic roots of these polynomials was 
done with Bini and Fiorentino's MPSolve package [4] which rigorously con- 
trols the output accuracy by dynamically altering the precision during com- 
putation. The computation of the largest real root was then cross-checked 
by using the bisection method with exact rational arithmetic in Maple. 

5 Roots Tending to Four 

To prove that a family of graphs X J 4 i ^(n) has real chromatic roots arbitrarily 
close to four, we will show that for any sufficiently small e > 0, there is an 
integer n such that for all n > no 

P(X AB (n)A-e) <0 and P(X AB (n), 4) > 0. 

We are mainly concerned with the situation when X A B(n) is planar, in which 
case the second inequality automatically holds. 
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The result of taking a matrix to a high power is essentially determined by 
its spectral properties. In this case, if < e < 2 — r, then for x = 4 — e, the 
matrix MD has four distinct positive eigenvalues 

Ai = 2 

A 2 = 2-5e+ 10e 2 /3 + O(e 3 ) 
A 3 = 2-8e + 26e 2 /3 + 0(e 3 ) 
A 4 = 

with corresponding (right) eigenvectors V\ = (1, — 1, — 1, 1) T , 

/ 3/2 + 35e/12 + 1103e 2 /216 + 0(e 3 ) \ 
1 + 5e/3 + 85e 2 /27 + 0(e 3 ) 
V2 ~ l + 5e/3 + 85e 2 /27 + 0(e 3 ) ' 

V - 1 / 

/ -e/3- 14e 2 /27 + 0(e 3 ) \ 
l/2 + e/6 + 4e 2 /27 + 0(e 3 ) 
V3 ~ l/2 + e/6 + 4e 2 /27 + 0(e 3 ) 

V 1 / 

and t>4 = (0, 1, —1, 0) T . The matrix MD is not symmetric, and so the eigen- 
vectors are not necessarily orthogonal. However, they do form an orthogonal 
set as long as we are willing to use a different inner product, namely 

(v, w) = v T Dw, 

where here, and subsequently, we identify the lxl matrix with the scalar it 
contains. 

5.1 Lemma. If {v±, . . . , t> 4 } are the eigenvectors of MD for some fixed x = 
4 — e where e G (0, 2 — r) then 

vJDvj = <^=^ i 7^ j. 

Proof The entries of D are all strictly positive and v i ^ so it follows that 
vj Dvi is a positive linear combination of squares with at least one non-zero 
term, and so vj Dvi > 0. Now as Vi and Vj are eigenvectors of MD we have 

{vjD)(MDv j ) = \ j vfDv j , 
(vjD)(MD Vl ) = XivjDvi. 
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The second expression is simply the transpose of the first, and so 

XjvjDvj = \vjDvj. 
If i ^ j, then Aj ^ Aj, and so it must be the case that vfDvj = 0. □ 

Let \\v\\ denote the norm of the vector v with respect to this inner prod- 
uct. As the eigenvalues {Aj} are distinct, the eigenvectors {v «} are linearly 
independent, so we can express Q(A) and Q(B) as linear combinations 

Q(A) = aivi + a 2 v 2 + a 3 v 3 + "4^4, 
Q{B) = (3 lVl + (3 2 v 2 + fcv 3 + fav 4 . 

Then the value of Q(A) T D(MD) n ~ l Q(B) at x = 4 - e is simply 

aiAAr'IMI 2 + a 2 p 2 \ n 2 - l \\v 2 \\ 2 + as^Arluair + onfaX^hA? ■ 

As n increases, this sum is dominated by the largest eigenvalue Aj for which 
both oti and $ are non-zero, and the eventual sign of this sum depends purely 
on whether the dominant ctj and $ have the same or opposite signs. 

We can now prove the result implied by the numerical computations of the 
previous section. 

5.2 Theorem. Graphs in the family X H ,w 4 (n) have real chromatic roots 
arbitrarily close to 4. 

Proof. As above, suppose that x = 4 — e with e G (0, 2 — r), that {v±, . . . , 1*4} 
are the eigenvectors for MD, and that Q(H) = Y and Q{W\) = Pi v i- 
Then substituting x = 4 — e into Q(H) T Dvi and Q(Wa) t Dvi for % — 1, 2 and 
expressing the results as series in e, we get 

ai||vi|| 2 = 0, 

AIMI 2 = o, 

a 2 |H| 2 = -50e + O(e 2 ), 
AdMI 2 = 5 + 20e/3 + O(e 2 ). 

Therefore the dominant eigenvalue for this family is A2 and for sufficiently 
small e > 0, the values a 2 and f3 2 have opposite signs. □ 
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6 Other end-graphs 



It is natural to consider which other end-graphs can be used in this construc- 
tion to produce real chromatic roots tending to 4. The graphs ^^(n) are 
planar if and only if both A and B are planar with a face of size 4 as the 
distinguished 4-cycle. In this situation, the coefficient of v\ in the expressions 
for Q(A) and Q(B) is always zero. 

6.1 Lemma. If A is a planar graph with a face o,iQ, 2 Q, 3 (i 4 of size 4 as distin- 
guished 4-cycle, then 

Q{A) T D Vl = 0. 

Proof. As v\ = (1, — 1, —1, 1) we need to prove that 

P 1 (A, x)/(x) 2 + P 4 (A, x)/{x) A = P 2 (A, x)/(x) 3 + P 3 (A, x)/(x) 3 . 

If x is an integer, then each term of the form Pi(A,x)/(x) s is the number of 
re-colourings of A of type % divided by the number of ways of choosing the s 
colours used on {ai, . . . , a 4 }. Equivalently, each term is the number of ways 
of extending a particular colouring of {a±, . . . , a 4 } to an re-colouring of the 
entire graph. Divide the colourings of A using colours from 1, 2, . . . , x into 
two classes as follows: p> is in Class 1 if 

[ip(a 1 ),ip(a 2 ),ip(a 3 ),ip(a 4 )\ = [1,2,1,2] or [1,2,3,4], 

and Class 2 if 

[ip(a 1 ),ip(a 2 ),ip(a 3 ),ip(a 4 )\ = [1,2,1,3] or [1,2,3,2]. 

We will now find a bijection between the two classes, and thus prove that 
the stated expression is true for all integers, and therefore for all x. 

Let ip be any colouring in either class, and consider whether or not there is 
a 2-4 path (i.e. a path whose vertices are all coloured 2 or 4) connecting a 2 
and 04. If there is a 2-4 path between a 2 and a 4 , then there is not a 1-3 path 
connecting a\ and a 3 because aia 2 a 3 a 4 is a face of a planar graph. Define 
a new colouring ip' by applying the following rules: if there is a 2-4 path 
from a 2 to a 4 , then exchange the colours 1 and 3 on the 1-3 component of A 
containing a 3 , and if there is no 2-4 path from a 2 to a 4 then exchange colours 
2 and 4 on the 2-4 component of A containing a 4 . It is straightforward to 
check that ip' is in Class 1 if and only <p is in Class 2, and vice versa. □ 
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We note that if A is planar, but the distinguished 4-cycle is not a face, then 
the conclusion of the above lemma may or may not hold. 

If C is a planar graph with a distinguished face of size 4, then say that it is 
positive if the expression for Q(C) T Dv2 at x = 4 — e has a positive leading 
term when expressed as a series in e, and negative otherwise. Then the family 
XA,B(n) has real chromatic roots arbitrarily close to 4 if one of A and B is 
positive, and the other negative. 

6.2 Lemma. If A is a planar graph with a face of size 4 as distin- 

guished 4-cycle, then the constant term in the expression for «2 11^2 1| 2 is 5/24 
times the number of 4-colourings of Aja\a^j a 2 a^ (i.e the graph obtained from 
A by identifying a\ with a 3 and a 2 with 0,4). 

Proof. When e = 0, the eigenvector v 2 = (3/2, 1, 1, -1) T = (5/2, 0, 0, 0) T + 
(-1, 1, 1, -1) T and so Q(A) T Dv 2 = 5/2 P^A, 4)/12 = 5/24 P X (A, 4). □ 

From this lemma, it follows that a necessary (but not sufficient) condition 
for a planar graph with distinguished 4-face to be negative is that the graph 
obtained by contracting the two diagonals of the face has no 4-colourings. 
Searching among the triangulations of a square with this property, negative 
graphs with as few as 10 vertices were found, with one example shown in 
Figure 6. 

6.1 Notes 

When Jacobsen, Salas and Sokal [12] examined the cylindrical triangular lat- 
tice without end-graphs, they obtained a 2 x 2 transfer matrix, with only 
A 2 and A 3 as eigenvalues. The reason for this is that rather than using 
colour partitions to index the rows and columns of the transfer matrix, they 
exploited planarity by using only "non-crossing non-nearest neighbour" par- 
titions. This is essentially equivalent to avoiding consideration of partitions 
of Type 1 which, as we have just shown, contain no additional information 
when A is planar. They gained another dimension by exploiting rotational 
symmetry, thus working directly in the subspace orthogonal to v±. Therefore 
the eigenvalues Ai = 2 and A4 = do not occur in their analysis. 

Rocek, Shrock & Tsai [14] obtained the same four eigenvalues for the transfer 
matrix as above, and also noticed that for certain pairs of end-graphs, the 
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Figure 6: A 10- vertex negative planar graph 



dominant eigenvalue is Ai = 2 (or rather, A4 = 2 in their notation) and that 
this caused the limiting curves of chromatic roots to change. In particular, 
the small closed curve passing through z = 4 and z ~ 3.481 becomes much 
larger, and passes through z = 4 and z = 1 + r ~ 2.618 instead (Fig. 4(b) of 
[14])- 
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